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Abstract 



• General formula for causal Green's function of linear differential operator of given 
degree in one variable, (d™ + Y^=o Pk(x)d^,), is given according to coefficient functions of 

| differential operator as a series of integrals. The solution also provides analytic formula 

^ for fundamental solutions of corresponding homogenous linear differential equation, (<9™ + 

Ew^W^lfW = as ser i es °f integrals. Furthermore, multiplicative property 
of causal Green's functions is shown and by which explicit formulas for causal Green's 
functions of some classes of decomposable linear differential operators are given. A method 
' to find Green's function of general linear differential operator of given degree in one 

• variable with arbitrary boundary condition according to coefficient functions of differential 
• \ operator is demonstrated. 

S ■ 

1 Casual Green's function for linear differential opera- 

cn : tors in one variable 

> 

Initial value problem for inhomogeneous linear differential equation of degree n in one variable, 

n-l 

{d x n + J2P^)d k x )y(x)=g(x), (1) 



can be converted to Volterra's integral equation of second kind [see [3]]. For initial condition 
d^y(a) — for i — 0, 1, .., n — 1, the corresponding Volterra's equation is given by; 

r x 

u(x) + / dzK(x, z) u(z) — g(x), 

J a 

where K(x,z) = (Y^ p k{x) ^=&]r ) and v( x ) = £ dz ^S)T u ( z )- For G 
L 2 [a, b], the condition ( dx dy \K(x, y)\ 2 ) < oo is sufficient condition for existence of unique 
solution in L 2 [a,b], given by iteration (e.g. see [2]). Clearly this conditions can be satisfied if 
Pi(x) ( i = 0, 1, n — 1) and g(x) functions are taken to be continuous on [a, b}. Therefore we 
can state the following theorem; 

Theorem 1 The Green's function for inhomogeneous linear differential equation (d™ + 
X]fc=o Pk{x) 9%) y(x) = g(x), where Pi{x) (i — 0, l,...,n— I) and g{x) are in C[a, b\, with 
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the boundary condition; d*y(a) = for i — 0, 1, — 1, is given by; 

G(x, y) = 0{x - y) ( ^~^ + f dz {X ~^ R{z, y)) , (2) 

where 

R(x, y) — h(x, y) + / ] I dz\ / dz% ■ ■ ■ I dz r -\ h(x, Z\) h(zi, Z/%) 
r =2 Jy 

••■h(zr-i,v), ( 3 ) 



n ~ 1 ■ \n-k-l 



and h(x, y) — — Pk{x)- — — -. The solution to inhomogeneous linear differential equa- 

tion (Op for x G [a, b] is then given by y{x) = P dzG(x, z) g(z). 



Proof. In order to prove ([2]) is Green's function of ([I]) it is enough to prove that for the 
two variables function; 

we have (d™ + J^Zl P k (x)d k x ) T(x, y) = and (d*T(x, y))\ x=y = for i = 0, 1, .., n - 2 and 
(<9™ _1 T(x, y))\ x=y = 1 [e.g. see [I]]. This can be easily done by noting; 

^' y)= ( (n-ll)! + J, dz \n-i-l)\ R(Z > y) ) i = Q ^- n ~ 1 ( 5 ) 

d^T(x,y) = R(x,y), 

and therefore, 

(d x n + £ T(s, y) = R(x, y) + J2 " ^ ^ 



(n-k- 1)! 

fc=0 fc=o v 7 

fc=0 v ' 

rx 

= R{x,y) — h(x,y) — / dz h(x, z) R(z,y) 
= R(x, y) - h(x, y) - (R(x, y) - h(x, y)) = 0. 

In the last line we used J y dz h(x, z) R(z,y) = (R(x,y) — h(x,y)^j, which comes from 
definition of R(x,y). By using (jSJ) we have (d£T(x,y))\ x=y = for i = 0,1,.., n — 2 and 
(dr 1 T(x,y))\ x=y = l. 

The Green's function for (CQ) with mentioned boundary condition is called causal solution 
which by method of variation of parameters is given by; 
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G{x , y) ^±3M^l )Hx - y) , (6 ) 



8=1 



where Ui(x),u 2 (x), ...,u n (x) are fundamental solutions of corresponding homogeneous dif- 
ferential equation; (d™ + Ylk=o Pk{x)d x ) Ui(x) = 0. W(y) is the Wronskian and Wi(y) is the 
Wronskian with its i th column in determinant is replace by (0, 0, .., 0, 1). Comparing this result 
with (j2J) we have the identity; 



vn-l px / _ z)n -l 



^ (n-l)! J y (n-l)\ V,y) V) 

For linear differential operator of first degree (n=l), like d x — P{x), the causal 
Green's function using [theorem [I] is equal to; (d x — P(x))^ 1 = 9(x — y) (\ + J2T=i Jy dzi 

■ ■ ■ f* n - 2 dz k . x J;*" 1 dz k P{z x ) ■ ■ ■ P{z k _ x )P{z k )) = 

0{X - V) (l + i EZl I;-'' Jy ' ' ' " ' ' P {z k -l) P {z k )) = 9{x - y)^ 

For linear differential operator of degree two in the form of; — P{x)), the causal Green's 
function by using [theorem [T] is given by T s (x, y)9(x — y) where; 

00 px rz k _ 2 rzk-i 

T s (x, y) = Ux -y) + y^A / dz 1 --- dz k -i / dz k 

k= i Jy Jy Jy 

(x - z l )P(z 1 )(z l - z 2 )P(z 2 )(z 2 - z 3 ) ■ ■ ■ (z k -i - z k )P(z k )(z k -y])\- (8) 



For example (d% — x) 1 = 9(x — y) (x — y) + 9{x — y) J x dz{x — z)z{z — y) + 9(x — 
y)r y dtp y dz ((x - t)t(t - z)z(z -„)) + ... = 9(x - y)((x - y) + - + M _ 

ll) + _ + ^72^ _ '^i^ + ^flo^ _ 504) + " ' )' wmcn is consistent with solution 



\ i „/ x /— Ai(x)Bi(w) + Ai(w)Bi(x)x , . , . 
(d? - x)' 1 = 9(x - y) ( t ,\ J , , — ; , / 7 r ) derived by @. 

1 x ' 1 A Ai( 2/ )Bi / (y) - Ai'(y)Bi(y) ; " 

It can be seen from (T5]) that if P%{x) ( z = 0, 1, n — 1) functions are smooth on [a, 6] then 
T(x,y), given by (j3J), is smooth function on [a,b] x [a, b], in which case we state the following 
theorem; 

Theorem 2 If Ti(x,y)9(x — y) and T 2 (x,y)9(x — y) are causal Green's functions for linear 
differential operators O x (x, d x ) = (<9™ + Yl=l p k{x)d x ) and 2 (x, d x ) = (d™ + YlT^o Qk(x)d^) 
respectively (assuming Pi(x) 's and qi(x) 's functions are in C°°[a, b}) then T 3 (x, y)9(x — y) where, 

T 3 (x,y)= dzT 2 (x,z)T 1 (z,y), (9) 
Jy 

is the causal Green's function for linear differential operator Os(x, d x ) = 0\{x, d x ).0 2 (x, d x ) 

Proof. By assumption; Oi(x, d x ) Ti(x,y) = and {d£Ti(x,y))\ x=y = for i = 0,l,..,n — 

2 and {d^T^y))^ = 1 also 2 (x, d x ) T 2 (x, y) = and {d*T 2 (x, y))\ x=y = for i = 
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0, 1, ..,m — 2 and (d™ 1 T 2 (x,y))\ x=y = 1, therefore we have; 

d*T 3 (x,y)= dztfifafazMT^y), 7 = 0,1,.., m-l (10) 



d™T 3 (x,y) = T 1 {x,y) + dz{d™{T 2 (x,z)))T x {z,y), (11) 



^T 3 (x,7;)=ar m (T 1 (x,7/)) + 9r m ( / cfe (^(T^,*)))^^)), 

Jy 

i — m + 1, .., m + n — 1. (12) 
Concentrating on the second term in ( fl2l) . we have for A; = 1,2, ..,77— 1; 



af( ^^(^(x^)))^^^) = {^^((a^-^c^^u 

Ti(x, y)) } + ( /; cfa (dr k (T 2 (x, z))) Ti(z, y)) 

= (ES Ei =0 ((5r +A - 1 - i+r T 2 (x, ,))i_ arn^, t/))} 

+(/;^(ar +fc (T 2 (x,,)))T 1 (,,7 / )) (is) 

From(riO|),(fTT]), (fT2|) and (H3J we have (<9jT 3 (z, y))\ x = y = for z = 0, 1, m+n-3, m+n-2 
and (9™ + " _1 r 3 (x, 7/))|x= y = 1- On the other hand; 

m-l »a; 

2 (x,9 x .)T 3 (x,7j) = (d x m + J2^(x)d x k ) dzT 2 {x,z)T x {z,y) 
= d x ( f dzd^ x T*(x,z)Ti{z,yj) 

/ x m—l 
dz(Y,Qk(x)dl:T 2 (x,z))T 1 {z,y) 
k=0 

rv 

= T 1 (x,y)+ dz0 2 (x,d x )T 2 (x,z)T 1 (z,y) 



X 



= Tt(x,y). 

Therefore O x (x, d x ).0 2 (x, d x )T 3 (x, y) = OtfadjTtfay) = 0. 

The following corollary comes as a consequence; 

Corollary 1 Causal Green's function for differential operator, 

0(x,d x ) = (d x -pi(x)) (d x -p 2 (x)) ■ ■ ■ (d x - p n (x)) , (14) 
where Pi(x) G C°°[a, b] (for i = 1, ...,n)is given by; 

G(x, y) = 6(x-y) f* dz x % dz 2 --- ^ ^-i ( A dU Pn(tn) 
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For example for differential operator 0(x, d x ) = d 2 + 3xd x + [2x 2 + 2), since (d x + x)(d x + 

2 

2x) = d 2 + 3xd x + (2x 2 + 2), by using result (JTSD one gets G(x,y) = J\ e y2 ~^ {Erf(^) - 
(Erf(^))} 9{x-y). 

Corollary 2 Causal Green's function for linear differential operator ; 

n 

0(x,d x ) = J2^kd x k , (16) 

where at GC and a n ^ 0, is given by; 

G(x, y) = / dzi / dz 2 - ■ ■ / dz n - 2 



Jy Jy Jy 



Zn-2 I \ 

dz i e l^ 1 (^c— ^i)+/32(^i— ^2) — h/3„(^„_i— y) ^) 



where f3±, {3 2 ■ • ■ /3 n are n complex roots of equation Q4X 1 = 0. 



i=0 



Proof. Differential operator 0(x,d x ) = Y^=o a i^x^ according to Fundamental theorem of 
algebra, can be written as, J2i=o a $x = a n(d x — Pi)(d x — /3 2 ) ■ ■ ■ (d x — j3 n ). Therefore by using 
( 1T5|) the result is proved. 

For example (<9 2 - u 2 )' 1 = 9(x - y) f y dz^*'^-^-^ = Bin M^) Q{x - y) and 
also (flg - zad 2 - w 2 ^ + zaa; 2 )" 1 = 9(x - y) dz^f* 1 d^ e M*-*0^fr-*a)+ta(*i»-v))) 

Q( x _ y) l e^-y^-e ia ^-y) _ sinh[a>(:r-y)] \ 

Lets consider a differential operator in form of 

0(x,d x ) = -dl+v{x). (18) 

By decomposing it into two firs degree differential operators; — (d 2 — v(x)) = —(d x —p(x))(d x — 
q(x)), we have consequently q(x) = —p(x) and p(x) 2 — d x p(x) = v(x). Therefore according to 
( IT5|) the causal the Green's function is given by; 

( - dl + ^a;))" 1 = -9(x - y) f dze^ -/WW dt '^ , (19) 



where p(x) is solution for first order nonlinear differential equation p (x) 2 — d x p (x) = v(x). 
This is just Riccati equation, thus the answers to p (x) 2 — d x p (x) = v(x) are given by solutions 

of homogenous differential equation (— d 2 + v(x))u\ j2 (x) = where p(x) = — {-^)- Inserting 

p(x) = — (^) into solution ffTQl . we have; ( — d 2 + v(x)) 1 = — 9(x — y)ui(x)ui{y) J* dz(-^-^). 

Considering the relation u 2 (z) = Ui(z) f dz — j—? (valid for homogenous differential equation 

(—d 2 + v(x))ui^(x) = 0) the solution (fl9l) becomes the standard solution, ( — d 2 + v(x)) = 



-9(x - y) \u 2 (x)ux(y) - ui(x)u 2 (y)j . 

Considering [theorem [2] one can introduce the following infinite non-abelian group of op- 
erators on a subspace of C°°[a,b]. We call it "Lalescu Group"; 

• Lalescu Group. The Group of differential operators of the form; (9^ + ^22=0 Pk( x )dx) °f 
all finite order, n>0, where Pk{x) G C°°[a,b} (for k=0,l,..,n-l) and their corresponding 
causal Green's functions G(x,y) = T(x,y)8(x — y) (given by on subspace ofC°°[a,b] 
consisting of functions which themselves and their derivatives to all orders are zero at 
x = a, creates non-abelian group with operators multiplication. 

Beside all differential operators 0(x,d x ) = [d x + Y^k=o Pk( x )d£) an d their causal 
Green's functions G(x,y) = T(x,y)9(x — y), the group also contains integro-differential op- 
erators and their inverses, coming from mixing these two groups of operators. For example 
Oi(x,d x ).T2(x,y)9(x — y), acting on <j){x) in the function space as 0±(x, d x )(f* dzT 2 (x, z)(f>(z)), 
and its inverse 2 (x,d x ).Ti(x,y)9(x — y) (where (x,d x ) = T 1 (x,y)9(x — y) and 2 l (x,d x ) = 
T 2 (x,y)9(x-y)). 



2 Green's function for linear differential operators in one 
variable with other boundary conditions 

In previous part the causal Green's function was derived according to coefficient functions of 
linear differential operator, here in this part it is shown that Green's function of general linear 
differential operator, for other boundary conditions on [a, b] (e.g. Sturm-Liouville problem), can 
also be derived according to coefficient functions of differential operator. First we note that ho- 
mogeneous linear differential equation of degree n in one variable, (d x + ^22=1 Pk( x )d x ) u(x) = 
0, for initial condition d*u(a) = q for i = 0,1,.., n — 1 can be converted to Volterra's 
integral equation of second kind (see [3]) as; (n(x) + f x dz K(x, z) fi(z)) = S(x), where 

K(x,z) = (rrJ P^) {j ^S^)> «(^0 = D(x)+IIdz^ff,(z), D(x) = EVo 2 ^ 

and S(x) — — Yl^o Ylk=o c i Pk ^%^k)\ — ■ Therefore we state the following theorem; 

Theorem 3 The Solution of homogeneous linear differential equation 
ipx + Ylk=o Pk{x) d x ) u(x) = 0, where Pi{x) (i = 0,1, ...,n — 1) are in C[a,b], with initial 
condition; d x u(a) = Ci for i — 0, 1, ..,n — 1, is given by; 

u(x) = D(x)+[ dzT(x,z)S(z), (20) 



where D(x) = y — , Six) = — } } Ci — and T(x,y) is given by Ol). 

fc=0 i=0 k=0 v ' 

It is easy to check that for boundary condition; d x u(b) = Ci for i = 0,l,..,n — 1 the 
answer for homogenous differential equation in C[a, b] is given by (J20l) in which a is replaced 
by b. One can derive solutions of homogenous differential equation, with specific boundary 
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conditions on either a or b, by using f[2"0"j) . By substituting solutions of homogenous differential 
equation with suitable boundary conditions (derived by ( 120]) ) in to the expressions for Green's 
functions given by method of variation of parameters, one can derive the Green's function 
according to coefficient functions of differential operator. For example the Green's function for 
Sturm-Lioville problem; 

(d 2 x -P(x))y(x)=g(x), B.C. y(a) = y(b) = 0, (21) 

by method of variation of parameters is given by; G(x,y) = (W(y))~ 1 (ui(x)u 2 (y) 0(y — x) 
+u 2 (x)ux(y) 6(x — y)), where U\(x) and u 2 (x) are answers of corresponding homogenous differ- 
ential equation, W(y) is the Wronskian and we have Ui(a) = u 2 (b) = 0. If we take B.C. Ui(a) = 

and d x Ui(a) = 1 then according to (TSUI) ; we have Ui(x) = (x — a) + f*dzT s (x,z)P(z)(z — a) = 
T a (x,a) (The function T s (x,y) is given by (jHJ)). On the other hand by taking B.C. 112(b) = 
and d x U2(b) = 1 we have u 2 (x) = T s (x,b). The Wronskian of U\(x) and u 2 (x) is constant and 
can be calculated easily at point W(y) = —T s (a,b) therefore; 

G(x, y) = (-T s (a, b))' 1 (t s (x, a)T s (y, b)9(y - x) + (x y)) . (22) 

3 Conclusion 

Causal Green's function for general linear differential operator in one variable was given by 
[theorem [T] as series of integrals. Multiplicative property of causal Green's functions is shown 
by [theorem [2] . For differential operators which are equal to multiplications of first order linear 
differential operators, explicit formula ffl5|) . was given for causal Green's functions. An infinite 
non-abelian group of operators on a subspace of C°° [a, b] is introduced. Analytic formula for 
fundamental solutions of homogenous linear differential equation in one variable was given 
via equation (JTJ) and for specific boundary condition via equation f )20|) as series of integrals. 
By using equation (F2UJ) and the method of variation of parameters a way to derive Green's 
functions with arbitrary boundary conditions in one variable according to coefficient functions 
of differential operators was given. 
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